Comparison of Measurements of Charge 
Transfer Inefficiencies in a CCD with 
High-Speed Column Parallel Readout 

Andre Sopczak 1 , Salim Aoulmit 2 , Khaled Bekhouche 1 , 
Chris Bowdery 1 , Craig Buttar 3 , Chris Damerell 4 , 
Dahmane Djendaoui 2 , Lakhdar Dehimi 2 , Rui Gao 6 , Tim Greenshaw 5 , 
Michal Koziel 1 , Dzmitry Maneuski 3 , Andrei Nomerotski 6 , 
Nouredine Sengouga 2 , Konstantin Stefanov 4 , Tuomo Tikkanen 5 , 
Tim Woolliscroft 5 , Steve Worm 4 , Zhige Zhang 6 

Lancaster University, UK 
2 Biskra University, Algeria 
3 Glasgow University, UK 
4 STFC Rutherford Appleton Laboratory, UK 
5 Liverpool University, UK 
6 Oxford University, UK 

Abstract 

Charge Coupled Devices (CCDs) have been successfully used in several high energy physics 
experiments over the past two decades. Their high spatial resolution and thin sensitive 
layers make them an excellent tool for studying short-lived particles. The Linear Collider 
Flavour Identification (LCFI) Collaboration has been developing Column-Parallel CCDs for 
the vertex detector of a future Linear Collider which can be read out many times faster than 
standard CCDs. The most recent studies are of devices designed to reduce both the CCD's 
intergate capacitance and the clock voltages necessary to drive it. A comparative study of 
measured Charge Transfer Inefficiency values between our previous and new results for a 
range of operating temperatures is presented. 
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Abstract — Charge Coupled Devices (CCDs) have been suc- 
cessfully used in several high energy physics experiments over 
the past two decades. Their high spatial resolution and thin 
sensitive layers make them an excellent tool for studying short- 
lived particles. The Linear Collider Flavour Identification (LCFI) 
Collaboration has been developing Column-Parallel CCDs for 
the vertex detector of a future Linear Collider which can be 
read out many times faster than standard CCDs. The most 
recent studies are of devices designed to reduce both the CCD's 
intergate capacitance and the clock voltages necessary to drive it. 
A comparative study of measured Charge Transfer Inefficiency 
values between our previous and new results for a range of 
operating temperatures is presented. 

Index Terms— LCFI, CPCCD, CCD, charge transfer ineffi- 
ciency, radiation damage 

I. Introduction 

The Nobel Prize-winning invention of an imaging semi- 
conductor circuit (the CCD sensor) [1], [2] has important 
applications for particle physics detectors. The study of ra- 
diation hardness is crucial for these applications [3]-[5], The 
LCFI collaboration has been developing and testing new CCD 
detectors for about 10 years [3]-[5]. Previous experimental 
results on CCD radiation hardness were reported for example 
in [6]-[10]. Several theoretical models have increased the 
understanding of radiation damage effects in CCDs [11]- 
[15]. Simulation and modeling of CCD radiation hardness 
effects for a CCD prototype with sequential readout was 
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reported at IEEE2005; comparing full TCAD simulations with 
analytic models was reported at IEEE2006; simulation and 
modeling of a CCD prototype with column parallel readout 
(CPCCD) was reported at IEEE2007 and in [14]. Experimental 
measurements using a method to determine the charge transfer 
inefficiency (CTI) were performed with a CPCCD prototype 
CPC-1 at a test stand at Liverpool University [17]. This work 
focuses on a new CPCCD prototype, CPC-T, at a test stand 
at Oxford University. The high radiation environment near 
the interaction point at a future Linear Collider damages the 
CCD material which leads to defects acting as electron traps 
in the silicon. The radiation level at a Linear Collider is 
estimated to be 5 x 10 11 e/cm 2 and 10 10 neutrons/cm 2 per 
year at the inner vertex detector layer (14 mm radius) [18], 
[19]. The mechanism of creating traps has been discussed in 
the literature [20]-[22]. These traps result in charge transfer 
inefficiency. 

The column parallel technology is in development to cope 
with the required readout rate. The CPC-T used is a 4-phase 
variant of the CPCCD technology capable of 50 MHz readout 
frequency. Experimental work at Liverpool University on an 
un-irradiated CPC-1 led to CTI values compatible with zero 
but with rather large uncertainties [17]. In this paper we 
demonstrate a method to determine the CTI value with an 
un-irradiated CPC-T aiming for small CTI uncertainties. 

II. Theory 

Soft X-ray photons (0.1 to 10 keV) interact with silicon 
atoms within the depleted layer. The depletion layer thickness 
is a parameter that determines the quantum efficiency at ener- 
gies above 4 keV [23]. The absorbed energy generates multiple 
e-h pairs. For a 5.9 keV X-ray source, one event (photon) 
generates a cloud of approximately 1620 electrons (Fig.[TJ[24]) 
contained within a diameter less than one micrometer [25]. 
The charge from a single X-ray photon, generated within the 
depletion region of a target pixel, is not transferred completely 
to the next pixel due to two main effects: the generation of 
thermal dark charge within the depletion region and the trap- 
ping of signal charge within the n-buried channel [26]. Since 
the buried channel is within the depletion layer, the important 
mechanisms are the capture of signal from the conduction band 
to the trap level and their subsequent emission back to the 
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Fig. 1. 55 Fe X-ray interacting with a CCD. A 5.9 keV photon generates an 
electron cloud of approximately 1620 e~ . 



conduction band [27]. Therefore, the X-ray event exhibits a 
'tail' of deferred charge. Also, the charge generated in the 
field-free region diffuses into neighboring pixels and adds to 
the 'tail' of deferred charge. The size and shape of this tail 
is a sensitive indicator of charge transfer inefficiency. X-ray 
stimulation is therefore extremely valuable in characterizing 
the CTI [24]. Many analyses have been made to simulate the 
effect of traps via the emission and capture processes [13], 
[26], [30]. The following simplified equations, based on earlier 
work by Shockley, Read and Hall [28], [29], have been used 
to analyse the CTI: 

dn t Tit + (N t - n t ) 

dt tv 
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where n t is the density of filled traps, Nt is the total density 
of traps, E t is the trap energy level below the bottom of the 
conduction band, r e is the emission time constant, t c is the 
capture time constant, a n is the trapping cross section, v t h is 
the thermal velocity of carriers, Nc is the effective density of 
states in the conduction band and n e is the density of electrons 
in the conduction band. For a detailed analytic model, the 
following parameters have been taken into account: 

• the order of magnitude of the emission and capture time 
constants compared to the shift time (time needed for a 
charge packet to move from one pixel to another). 

• the shape of the electrostatic potential, which can be 
assumed to be placed in the middle of the well. 

• the level of the signal charge (density of free electrons) 
within the potential well in comparison to the total density 
of traps. 

III. Test Stand for CCD Operation 

A test stand has been set up with readout electronics and a 
freezer unit as shown in Fig. |2] The temperature range of the 
freezer is from room temperature down to about —60 °C. Fine 
control of the CPC-T temperature is done using a CAL9500P 



controller (the temperature is kept constant within 0.1 °C). 
A flux of boiled nitrogen is introduced into the motherboard 
box to purge water vapour. The CPC-T chips come in 2 main 
variants: inherent 4-phase CCD driven as 2-phase CCD, and 
'pedestal' 2-phase CCD with 2 additional DC-biased gates. 
The former was used for this measurement. The first and 
second gates of each pixel, P1A and P1B, are driven by 
Phase 1, and PI A is offset by the DC voltage OPV (Offset and 
Pedestal Voltage), as shown in Fig. [3] The CPC-T has 500 x 10 
pixels with a pixel size of 20 x 20 /im 2 . Initial measurements 
have been performed on an un-irradiated device in standalone 
mode, where the signals from four columns of the CCD were 
amplified and connected to external 14-bit ADCs. A 55 Fe 
source emitting 5.9 keV X-rays was attached to a holder at a 
distance of 1 cm from the CCD to provide the signal charge. 
The schematic diagram in Fig.|4]illustrates the electronics used 
to drive and read out the CPC-T. The apparatus is controlled 
by a Lab View program through interface modules. The BVM2 
sequencer receives the master clock of 1 MHz from the 
function generator to provide four signals, two for ADC and 
two to trigger the generators which produces a CCD clock 
and reset gate signals. The 2-resets configuration, when one 
reset is applied before reading the first pixel of the CCD and 
one after reading the last pixel, is used in this measurement. 
This configuration leads to low noise since the reset noise is 
absent. The occupancy is about 1% for the integration time 
of 100 ms given the strength of the X-ray source and the 
experimental layout. The number of frames (complete readout 
of the CCD) is not kept constant to study the effect on the 
statistical uncertainty. Our method is based on the typical 
methods used for serially read out CCDs, where the CTI is 
determined by fitting a line to the readout charge signal as a 
function of the pixel number. A linear function can be expected 
when the CTI is small. 
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Fig. 2. Picture of the CPC-T readout. The CPC-T mother board is inside 
the freezer. 
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Fig. 3. Clocking scheme for the 4-phase variant. This variant is driven as a 
two-phase CCD with additional DC voltage (OPV) to the voltage of the first 
gate of a pixel. 
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Fig. 4. Schematic diagram of the CPC-T readout. 



IV. Signal Measurement and CTI Determination 
Method 

The fast ADCs convert the signal charge after amplification 
with a wideband preamplifier. The four columns are read out 
in 4 channels by ADCs. The signal charges of 500 pixels 
were acquired in 1000 and 10000 frames per measurement for 
each temperature. Using 10000 frames leads to a sufficient 
statistical precision (around few times 10 -6 for the CTI 
measurement). The collected data have been analysed using 
MATLAB [16]. First, we begin by applying correlated double 



sampling, where the difference between the signals of two 
consecutive pixels is taken to be the signal charge collected 
by the latter pixel. This reduces some of the noise components 
(e.g. 1/f, kTC, white noise, etc.) in the CCD signal. As an 
example, Fig. [5] shows the pulse-height distribution of ADC 
codes for the last 10 adjacent pixels in column 2. The charge 
transfer inefficiency (CTI) in one pixel is defined as the ratio 
of signal lost during transfer (captured by traps) to the initial 
signal charge. The CTI is calculated following these steps: 
creation of a histogram with ADC codes of 10 pixels in 
a column. These pixels have nearly the same baseline. The 
histogram creation is repeated 50 times to cover all 500 pixels 
in a column. Fits with Gaussian functions are made to the noise 
and X-ray peaks. We use xo — na and xq + na as limits for 
noise peak and X-ray peak (xo and a are the centroid and the 
standard deviation respectively of the two Gaussian functions 
resulting from fitting the noise and X-ray peaks) to determine 
noise centroid and X-ray centroid for each pixel. The factor n 
is chosen between 1 and 3 depending on the amount of charge 
sharing. The X-ray signal for each pixel is the difference 
between the X-ray peak centroid and the noise peak centroid. 
Figure [6] shows the distribution of the averages as a function 
of pixel number. The distribution is fitted with the first-order 
polynomial function Po + Pij, where Pq corresponds to the 
charge at the first pixel, Pi is the slope and j is the pixel 
number. The CTI is determined using CTI = —Pi/Pq. 
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Readout frequency=1 MHz, 
T=-27.5 °C, OPV=-0.81 V, CLK voltage=4 V. 
integration time=100 ms, 10000 frames. 
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Fig. 5. Distribution of ADC codes for channel (column) two. Three regions 
are observed: (I) the high peak region which represents the noise, (II) the 
region separating the two peaks which represents the charge sharing between 
pixels, and (III) the X-ray peak region which represents the fully collected 
charge in a single pixel. The noise and X-ray peaks are fitted by a Gaussian 
function to determine the centroids. The X-ray signal is the difference between 
the two centroids. 
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Readout frequency=1 MHz, 
T=-27.5 °C, OPV=-0.81 V, CLK voltage=4 V. 
integration time=100 ms, 10000 frames. 
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Fig. 6. Linear fit to average ADC codes. The X-ray centroid of each pixel is 
calculated by averaging ADC codes within the interval xq — no and xq + no- 
where Kn<3. 
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Readout frequency=1 MHz 
integration time=100 ms, 
OPV=-0.81 V, 
CLK voltage=4 V. 
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Fig. 7. CTI as a function of temperature for different numbers 
The error bars have been significantly reduced by increasing the 
frames. 
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V. CTI Results Pre-Irradiation 

Figure [7] shows the CTI values as a function of temperature 
for an un-irradiated CPC-T for 1000 and 10000 frames. The 
CTI has been calculated using a linear fit of average ADC 
codes versus pixel number. Figure [8] shows a comparison with 
the CPC-1 measurement [17] taken at the test stand in Liver- 
pool. Uncertainties have been reduced mostly by increasing the 
number of frames. For this CCD with 500 pixels per column a 
CTI value of 10~ 5 means that only 0.5% of the signal charge 
is lost, which is acceptable in normal operation. The apparent 
trend of the CTI at high temperatures in the operating range 
used is probably due to the contribution of two effects. First, 
there is the effect of thermal carrier generation (dark current) 
which is highly temperature-dependent. The dark current, non- 
uniform by nature, can have a large effect on the signal charge 
transfer for high temperatures, long integration time and large 
number of pixels in the column [31]. Second, there is the 
possibility of the presence of low trap density that could have 
been created during the long duration (around two years) of 
exposure to a soft X-ray source while studying the device. 
This significant positive value of CTI before irradiating the 
CCD was observed experimentally and modeled by a simple 
analytic model by including one trap level [4]. Using our 
analytic model [15], the CTI is expressed as 

CTI = 2^ [1 - exp(-t(- + -))] x 

U T S (l-CXp(-^)) * 

-exp( — -)]. 

We have fitted the CTI curve including two deep traps as 
shown in Fig. [9] Both are electron traps at 0.37 eV and 0.44 eV 
below the bottom of the conduction band and having a trapping 
cross-section a n = 3 x 10 -15 cm 2 . We have considered that 



there is no interaction between the two traps, so they interact 
independently with the signal charge. Therefore, the total CTI 
is the sum of CTIs resulting from the effect of the traps. The 
fit is in good agreement with the data and shows that the 
0.44 eV trap is the dominant one in this range of temperature 
as its density (TV^ = 5.22 x 10 10 cm -3 ) is much larger than 
that of the 0.37 eV trap (N a = 2.63 x 10 9 cm" 3 ). 

The accuracy of the CTI calculation can be improved by 
positioning the 55 Fe source so that it irradiates uniformly 
the CCD, carefully choosing the gain to use the ADC in its 
maximum dynamic range and acquiring data in a large number 
of frames. The non-uniformity of the X-ray source coverage 
has an effect on the CTI determination. This is well understood 
and reproduces the estimate based on a geometry where the 
source is placed 1 cm away from the CPC-T. Figure [10] 
shows the non-uniformity of the X-ray source coverage. The 
figure contains two curves, the measured and the estimated 
X-ray distributions. The measured distribution is determined 
by counting all ADC codes above the noise threshold Xq + 3a. 
The estimated distribution is determined using the following 
formula for the given geometry: 



f(n) = 



(h/l p ) 



(h/l p ) 2 + (n - n )' 



3/2 



where / is the distribution of X-rays upon the CPC-T, n is 
the pixel number, h is the distance between source and CPC- 
T, l p is the length of one pixel and no is the pixel number 
corresponding to the vertex position. In order to avoid the 
effect of the non-uniformity, the first and last 50 pixels are 
excluded from the fit to the averages for the CTI determination 
(Fig.0. 
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Fig. 8. Comparison of the Oxford results with 10000 frames with the 
Liverpool results [17] where the number of frames was 5000 and a fraction 
of data was lost because of a sampling inefficiency. 
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Fig. 10. X-ray distribution upon CPC-T. Measured distribution is determined 
by counting all ADC codes above the noise threshold which represent the 
energy deposited by X-rays. Simple formula reflecting the geometry is used 
to estimate the distribution. 
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Fig. 9. Non-linear fit of the measured CTI using our analytic model [15]. The 
model includes two acceptor traps, 0.37 and 0.44 eV below the conduction 
band. A trapping cross-section a n = 3 X 10" 15 cm 2 is used for both traps. 

VI. Conclusions and Outlook 

An un-irradiated CPC-T was operated in a range of tem- 
peratures from —15 °C to —60 °C (freezer cooling) with 
different numbers of frames, 1000 and 10000. The CTI is 
analysed at different operating temperatures. A clear X-ray 
signal is extracted by calculating the difference between the 
noise centroid (baseline) and the X-ray centroid after well 
determining them. The statistical uncertainties have been re- 
duced compared to a previous work with CPC-1 [17]. The 
reduced uncertainties are due to the improved method, as we 
do not expect a different behavior between CPC-1 and CPC-T 
regarding the CTI measurements as both devices have a very 
similar geometric structure. The measurement of significant 



non-zero CTI values are indicators of trapping and thermally 
generated electrons. The former is dominant in this range of 
low temperatures. We expect that the CTI will increase after 
irradiation as the trap density increases. The non-uniformity of 
the X-ray radiation has to be taken in account when measuring 
the CTI of an irradiated CPC-T in the future. 
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